Single hole dynamics in the t-J model on two- and three-leg ladders 
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The dynamics of a single hole in the t-J model on two- (2LL) and three- (3LL) leg ladders is 
studied using a recently developed quantum Monte Carlo algorithm. For the 2LL it is shown that 
in addition to the most pronounced features of the spectral function, well described by the limit of 
strong coupling along the rungs, a clear shadow band appears in the antibonding channel. Moreover, 
both the bonding band and its shadow have a finite quasiparticle (QP) weight in the thermodynamic 
limit. For strong coupling along the rungs of the 3LL, the low-energy spectrum in the antisymmetric 
channel is similar to a one-dimensional chain, whereas in the two symmetric channels it resembles 
the 2LL. The QP weight vanishes in the antisymmetric channel, but is finite in the symmetric one. 
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The n-leg ladder systems serve as a bridge from one- 
dimensional chains to two-dimensional layered cuprates, 
offering thus a rich playground for studying the interplay 
of charge and|-s»in degrees of freedom in strongly corre- 
lated systemsEra. In particular, two- (2LL) and three- 
(3LL) leg ladder materials became recently n available, 
making experiments possible in such systemstj. At zero 
doping, the spin-excitation spectrum is gapped in the 
case of even-leg ladders, whereas it is gapless in the case 
of an odinumber of legs, as found by experiments on 2LL 
and 3LLotl and by quantum Mcate Carlo (QMC) simu- 
lations of up to five-leg ladders! A discussion on the 
impact of this odd-even effect on the charge excitation 
spectrum as obtained by angle resolved pJjiotoemission 
(ARPES) experiments has recently startedU. 

Charge carriers on cuprates can be well described by 
the t-J modeltl For the 2LL a considerable amour 
work has been done by exact diagonalizations (ED)E 
density-matrix renormalizatian group (DMRG)cS, and 
strong coupling expansionsEilE3, leading to the following 
picture: When the coupling along the rungs J± is much 
larger than along the chains (Jn), the bonding and the 
antibonding band are split by 2ij_. The bonding band 
shows a sharp peak at the lower edge of the spectrum. 
Approaching the isotropic case, the energy gap between 
the bonding and antibonding band is reduced, and the 
bands become nearly flat on half of the Brillouin zone. 
Results of exact diagonalizationsEa indicate, that spin and 
charge excitations are coupled to each other, leading to 
the existence of quasiparticles in these systems. Like 
for the 2LL, the strong-coupling limit of the 3LL is a 
promising starting point. In this limit, the low-energy 
excitations have odd parity, and can be described by an 
effective one-dimensional (ID) t- J-modelli 3 ], which forms 
a Luttinger liquid in a wide parameter range of J/t. Re- 
sults from ED of three-leg ladders suggest that the two 
even channels have the same physics as 2LL|j up to the 
isotropic case. However, due to the small systems con- 
sidered, the gap between the even and the odd channel 



is of the same order as the finite-size spin-gap. Whether 
a particular channel will evolve upon doping towards 
the Luttinger- Tomonaga or the Luther-Emery universal- 
ity class is signaled by the absence or not of a quasiparti- 
cle (QP) weight in the single hole case. This can only be 
determined by performing a finite-size scaling. Whereas 
a Mott-insulator with branch cuts in the spectral func- 
tion (i.e. without QP) evolves towards a Luttinger-liquid 
upon doping, the Luttinger-liquid parameter reaches the 
universal value K p p. 1 in the limit of zero doping for a 
spin-gapped systemEJ. Then, using a previous analysis of 
the spectral function for Luther-Emery systemal3 a finite 
QP weight in the single hole case results. 

In this paper we study single hole dynamics in the 
t-J model on 2LL and 3LL using a recently developed 
algorithm, which has-been successfully applied to 



oneta and two dimensionsEZl. The spin.-dynaniics is ef- 
ficiently simulated by a loop-algorithmlia and for each 
spin configuration, the one-particle Green's function is 
calculated exactly. For details see Ref. 17. The QP 
weight and dispersion are determined from the asymp- 
totics in imaginary time of the Green's function and the 
full spectral function is obtained with the Maximum En- 
tropy methocO. 

The t-J model is given by 



Jij (Si ■ Sj /n-iUj) , 



<i.3> 



(i) 



where cj are electron operators restricted to the 
Hilbert space with no double occupancy, Si — 



( 1 / 2 )Ea,/ 3 4a a;Q ./3 c i,/9' = EjCjc^ and the sum 
runs over nearest neighbors only. For i,j on a rung we 
set J^ — J±,tij = t±, and Jy = J\\,tij = t\\ if i,j are 
on the same leg ( t = tii is the unit of energy in the 
following) . Along the legs of the ladder we use periodic 
boundary conditions, whereas the boundaries are open 
along the rungs. The reference energy is the Heisenberg 
groundstate. 
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We start our discussion with the 2LL in the limit 
J± 3> J||. At half filling, this leads to the formation 
of singlets on each-jungE]. Perturbative treatments^! or 
series expansionsLHI around this limit show that, like in 
the free case, the hole has a cosine-like dispersion with a 
splitting of the bonding and the antibonding bands by t± . 
The only difference is a rescaling of the effective hopping 
to — 1| | / 2 and a shift of the energies by J±. 
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FIG. 1. Quasiparticle dispersion and spectral function for 
the bonding (a) and the antibonding (b) band for the 2LL 
with J±/t\\ = 1.6, J\\/t\\ = 0.4, and t x /t\\ = 2. Further 
details are discussed in the text. 

Figure [j] presents our results for a 2LL with 2 x 32 
sites, where J\\/t\\ = 0.4, Jj_/i|| = 1.6, and t±/t^ = 2. 
For clarity, the quasiparticle peaks were cut off at a given 
intensity^. As expected, the major feature in the anti- 
bonding channel is 4i = 2t± above the bonding band. 
Both bands can be fitted by (0.75 cos — 0A7)t for the 
bonding and (0.75 cos fen — 0.47 + A)t in the antibonding 
band (full lines in Fig. |l|). However, we find addition- 
ally a clear evidence of a shadow band in the antibond- 
ing channel. This is corroborated by considering apart 
from the spectral function, the QP dispersion directly ob- 
tained from the imaginary time Green's function (x with 
errorbars), and by superposing the fitted bonding band 
with an [ energ y shift of 0.64i, which is approximately the 
spin gapE3'Eil, and a momentum shift of 7r [dashed line in 
Fig. 0(b)], as expected for a shadow resulting from the 
coupling of the bonding band to the lowest spin excita- 
tions centered at (7r, tt) and with an energy given by the 
spin gap. The possibility of a shadow band for 2LL was 
first proposed on the basis of EDlI, but due to the small 
system sizes, no quantitative assignment was possible. In 
principle, also a shadow of the antibonding band should 
be expected in the bonding channel. However, since the 
spectral weight in the antibonding channel is split be- 
tween the original band and the shadow of the bonding 
band, resulting in rather broad features in comparison to 
the bonding band, less well defined shadows should be ex- 
pected in this case. Only a weak structure is observed in 
Fig. 0(a) at around the energy of the antibonding band, 
but it is difficult to assign a dispersion to it. 

Figure H presents the isotropic case Jm = Jj_ = tii = 
for the 2LL (2 x 48 sites). rAs-before, expansions around 
the strong coupling limitll^'ta describe accurately the 



bonding band obtained from our simulations. The high- 
energy excitations show no clear dispersion, and are 
about 5t higher in energy in both bands. Again, the 
antibonding band can be interpreted as a shadow of the 
bonding band. However, the energy difference between 
the bonding band and its shadow (A = 0.17t ± 0.02iVis 
now smaller than the value of the spin gap (~ J/2)E3, 
which is only an upper bound that is reached when the 
bonding and antibonding bands are well separated as in 
the strong coupling case. 




FIG. 2. Quasiparticle dispersion for the bonding (a) and 
antibonding (b) band for the isotropic two-leg ladder with 
J/t — 1. Further details are discussed in the text. 

The lowest states for the hole in both the bonding and 
its shadow are now closer to k\\ — n/2 with a band- 
width ~ 0.5t, ip_.rough agreement with ARPES results 
on Sri4Cu2404itj in spite of the somewhat large value of 
J. The bands observed in ARPES are symmetric around 
fc|| = 7r/2, a feature that can be reproduced by super- 
posing the bonding band and its shadow in Fig. ||. The 
very flat portions around fc|| = ir (k\\ — 0) for the bond- 
ing (shadow) are reminiscent of the flat dispersion around 
(7r, 0) or (0, 7r) in two dimensions, as one approaches these 
points from (0, 0) or (7r, ir) respectively, as already sug- 
gested based on a reduced basis approximations. 

We consider next the 3LL and start with the strong 
coupling limit, where the low energy behavior was shown 
to correspond to an effective piie-dimensionaLt- J-chain 
in the antisymmetric channelta, whereas EDB indicate 
that the two symmetric channels have a finite energy 
gap to the antisymmetric one. The observed gap is how- 
ever quite small ~ 0.3J = 0.15i, and therefore ^t^tf-i 
size effects cannot be excluded. Results from DMRGIIJE^I 
show, that at low doping, the holes are at the two outer 
legs of the 3LL confirming the picture derived from ED. 
Figure || shows the spectral function for a 3LL with 
3 x 32 sites and the same parameters as for the 2LL. 
As expected, the lowest state belongs to the antisym- 
metric channel. For Jj_ 3> J\\ , the effective parame- 
ters of the corresponding t-J chain are given by t e ff = 

3t|| (j± + y/J'i + 8t^j /S^J'i + St 2 ± and J eff = J||E1. 
Using these-narameters, and a model with free spinons 
and holonsE3ll3, a reasonably good description of the 
lower edge of the antisymmetric band is obtained (full 
line in Fig. |||(a) corresponds to the minimal energy given 
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by the convolution of spinons and holons). Even more in- bands. They are well separated from the antisymmetric 
teresting are the results obtained for the two symmetric one in energy, and can be very well fitted by the bonding 




Ujt u/t Uj/t 

FIG. 3. Spectral function for the 3LL with J±/t| = 1.6, J\\/t\\ = 0.4, and t±/t^ — 2 in the antisymmetric (a), the symmetric 
bonding (b) and the antibonding (c) band. Further details are discussed in the text. 




FIG. 4. Spectral function for the 3LL with J/t = 1. Further details are discussed in the text. 



band dispersion for the isotropic 2LL from perturba- 
tion theoryliil, with J/t — 0.4, and its shadow [full lines 
in Fig. ||(b) and (c)]. The two bands with even par- 



ity are connected. 
which are gapless& 



spin excitations with even parity, 
such that in this case no shift be- 
tween them appears. The description of the 3LL system 
at the isotropic point (J = t,t e ff = 0.5,3 x 32 sites, Fig. 
|j) on the other hand, is less satisfactory when the strong 
coupling limit is used. All three bands have their min- 
ima around 37r/4, and contrary to EDO, we do not find 
a significant energy gap between all three bands. The 
shape of the lower edge of the antisymmetric channel - 
which is determined by the holon dispersion relation for 
this value of J e / / / t e f / = 2 - differs considerably from 
that of the Bethe-Ansatz_holon dispersion at the super- 
symmetric point J/t — 2a. 



address the QP weight Z(k) = 
2 which can be obtained from the 



We now 

imaginary time Green's function as the weight of the 
exponential with the slowest decay at large rt. As can 
be seen in Fig. ||, the QP weight is finite in the ther- 
modynamic limit for both the bonding band and its 
shadow in the 2LL, and large compared to the 2D t-J 



modelEl (this result is valid for all fc-points where the 
data is accurate enough to extract the QP weight). On 
the contrary, finite-size scaling of the QP weight in the 
antisymmetric channel of the 3LL leads to a vanishing 
QP weight in the thermodynamic limit, as expected for 
a Luttinger liquid. The scaling of Z(k) can be fitted 
by a * L b with a = 0.464 ± 0.014,0.564 ± 0.015 and 
b= -0.335 ±0.015, -0.306 ±0.016 for J/t = 1,2 respec- 
tively. Squares in Fig. || (c) and (d) correspond to ED 
results that confirm the QMC ones. The values obtained 
for the exponent b depart from what would be obtained 
by naively using the strong coupling result, that for J = t 
leads to t e ff = 0.5 (i.e. the supexsymmetric case) where 
the exponent should be b = — 0.5E1 In contrast to the re- 
sults for the antisymmetric channel, the QP weight stays 
finite for the symmetric channel in the thermodynamic 
limit [Fig. [51(b)], although smaller than in the 2D caseO. 
Therefore, our results demonstrate that the antisymmet- 
ric channel in the 3LL will evolve upon doping towards 
the Luttinger- Tomonaga universality class, whereas the 
symmetric channel belongs to the Luther-Emery one, as 
suggested previouslytj. 
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FIG. 5. Finite-size scaling for 
fcn = 7r/2 for the bonding band 
(k± = 7r) in the 2LL at J/t = 1.2 
nel of the 3LL (b) at J/t — 1 and 
bonding band (full symbols) and 
ing band (open symbols)], and in 
the 3LL at fe|| = tt/2 for J/t = 
The squares are results from ED 



the quasiparticle weight at 
(k± — 0) and its shadow 
(a), in the symmetric chan- 
J/t = 2 [fc|| = 3tt/4 for the 
fc| = it /A for the antibond- 
the antisymmetric band of 
1 (c), and for J/t = 2 (d). 



Summarizing, wc have studied the spectral properties 
of a single hole in the t-J model on 2LL (up to 2 x 96 
sites) and 3LL (up to 3 x 64 sites) using a recently de- 
veloped QMC algorithm. It is shown, that expansions 
around the limit J± 3> J\\ describe accurately the bond- 
ing band and the high energy portion of the antibonding 
band in the 2LL, even in the isotropic case. However, 
such an expansion misses the low-energy portion of the 
antisymmetric channel, that corresponds to a shadow of 
the bonding band. For the 3LL strong coupling expan- 
sions (Jj_ 3> Jii) for the antisymmetric channel and the 
2LL give a good description of the antisymmetric and 
symmetric channels respectively, as long as J± is appre- 
ciably larger than Jy. The dispersions in the isotropic 
3LL seem to be less well described by a strong coupling 
expansion. The QP weight extrapolates in the thermo- 
dynamic limit to a large finite value for the bonding band 
and its shadow in the 2LL. For the 3LL finite-size scal- 
ing leads to a vanishing QP weight for the antisymmet- 
ric channel but a finite one for the symmetric channel, 
demonstrating that the first one corresponds to the Lut- 
tingcr universality class, whereas the latter corresponds 
to the Luther-Emery one. 
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